P153  24?  NUMERICAL  METHODS  FOR  STIFF  ORDINARY  AND  ELLIPTIC 
PARTIAL  DIFFENENTIAL  EG  (U)  IBM  THOMAS  J  WATSON 
RESEARCH  CENTER  YORKTOUN  HEIGHTS  NY  L  WERNER  ET  RL. 
LASSIFIED  FEB  85  AFOSR-TR-85-0336  F49620-82-C-0005  F/G  12/1 


Ilbu 


AFOSR-TR 


~  *  0  3  3  6 


NUMERICAL  METHODS  FOR  STIFF  ORDINARY  AND  ELLIPTIC  PARTIAL 
DIFFERENTIAL  EQUATIONS 


Final  Report 


I.  SUMMARY  OF  PROGRESS 


DTIC 

f^ELECTE 
APR  30  885 


$ 


B 


1.  Wave-form  (WF)  relaxation  methods  for  large-scale  circuit  analysis 


Work  has  been  completed  or  has  progressed  on  items  2b,  2c,  and  2e  of 
Section  IV  of  the  Contract  Proposal.  More  specifically:  The  convergence 
of  the  discretized  version  of  the  VR  algorithm  was  shown  [1]  wider  suit¬ 
able  assumptions  on  the  stability  of  the  multistep  methods  employed  and 
on  the  strength  of  the  feedback  (item  2b).  Also,  the  effects  of  various 
decompositions  of  the  large-scale  circuit  were  discussed  [1,2]  (item  2e). 
In  particular,  in  [1],  a  new  decomposition  which  uses  only  neighboring 
circuits  was  shown  to  be  effective  for  a  large  class  of  digital  circuits. 
This  decomposition,  called  t -scheduling,  dramatically  decreases  the 
storage  requirements  of  the  VR  method.  Moreover,  since  different  sub¬ 
circuits  are  usually  integrated  using  different  time  steps,  the  effect 
of  interpolation  on  the  convergence  of  the  VR  method  was  studied  (item 
2c).  Various  sufficient  conditions  for  the  VR  convergence  were  derived 
which  combine  accuracy  requirements  with  some  sort  of  stability  condi¬ 
tion,  e.g.  that  the  interpolation  operator  be  of  positive  type,  or  be 
"passive",  i.e.  act  boundedly  on  exponentially  increasing  sequences. 
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The  above  theoretical  results  have  been  used  in  a  program,  for  hierar¬ 
chical  analysis  of  large  scale  MOSFET  circuits,  which  implements  an  ef¬ 
ficient  scheduling  (partitioning)  algorithm  (item  3e)  and  uses  the  charge 
as  a  state  variable,  leading  to  better  numerical  stability.  The  program 
is  described  in  reference  {3]. 


2.  Fast  direct  methods  for  elliptic  partial  differential  equations 

Work  has  been  completed  or  has  progressed  on  items  3a-3d  and  3f  of  Section 
IV  of  the  Contract  Proposal.  More  specifically:  A  one-parameter  family 
of  factored  complex  discretizations  of  the  Laplace  operator  was  derived 
[A].  These  discretizations  are  second  order  accurate  and  serve  as  a  basis 
for  fast  direct  or  iterative  Poisson  solvers  on  general  two-dimensional 
regions  (item  3a).  The  stability  of  these  discretizations  with  respect 
to  the  propagation  of  rounding  errors  was  analyzed  [A].  A  variant  of  the 
marching  method,  based  on  these  factored  discretizations,  was  proposed 
[A]  which  is  much  more  stable  than  the  conventional  marching  method  and 
is  thus  applicable  to  grids  with  large  numbers  of  discretizations  steps 
in  each  direction,  without  resorting  to  domain  decomposition  or  multiple 
shooting.  The  gain  in  stability  is  achieved  by  solving  initial  value 
problems  for  two  first-order  difference  equations,  rather  than  one 
second-order  equation.  Furthermore,  the  stability  and  accuracy  of  the 
method  can  be  controlled  by  the  choice  of  the  above  mentioned  parameter 
(item  3b).  The  new  version  of  the  marching  method  was  implemented  and 
successfully  tested  [A]  on  a  square  region  as  a  function  of  the  number 

of  discretization  steps  (item  3c),  and  on  several  examples  of  general 

AIR  fOVfS  . -  ?  "  T  r  *1"  T  FT  T  '  ***'*' 

HOT  I". 


-  • 

D  t  :i  t  r  : 

MS iTHi  :  7- 
Chief,  T c c r. • 


•  i  i  i  v  ; 


f  on 


regions  with  piecewise  rectilinear  boundaries,  including  multiply  con¬ 
nected  ones  (item  3d).  The  numerical  results  confined  the  theoretical 
accuracy  and  stability  predictions.  The  aarching  method  was  extended  to 
the  case  of  general  nixed  linear  boundary  conditions  (itea  3f),  and  the 
stability  analysis  was  generalized  to  the  use  of  rectangular  (rather  than 
square)  discretization  neshes. 

II.  PROPOSED  RESEARCH 

During  the  second  year  of  the  contract  period,  we  propose  to: 

la)  Investigate  the  effects  of  stiffness  on  the  convergence  of  the 
VR  aethod,  for  both  continuous  systems  and  their 
discretizations. 

lb)  Devise  algorithms  for  detecting  strong  feedbacks  in  large  cir¬ 
cuits  and  strategies  for  "windowing"  VR  in  the  presence  of  such 
feedbacks,  with  the  object  of  optimizing  the  over-all  effi¬ 
ciency  of  the  method. 

lc)  Carry  out  numerical  testing  of  the  VR  method. 

2)  Investigate  the  short-channel  effects  in  MOS  transistors  by 

both  analytical  asymptotic  methods  and  by  numerical  integration 
of  the  partial  differential  equations  governing  and  typical 
semi-conductor  devices. 
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Study  the  us«  of  th*  pr* -conditioned  conjugate  gradient  (PCG) 
•at hod,  in  conjunction  with  th*  factored  consistent  complex 
discretisations  of  th*  Laplacian  Motioned  above,  for  th*  it¬ 
erative  solution  of  Poisson's  equation  on  general  two- 
dimensional  regions. 

Apply  the  PCG  aethod  to  approximate  factored  discretisations 
of  th*  Laplacian,  with  a  view  toward  stabilising  the  computa¬ 
tion  and  deriving  iterative  algorithms  which  are  faster  than 
th*  direct  algorithm  dascribed  above. 
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